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Abstract 

Using an averaged generating function for coloured hard-dimers, some ran- 
dom variables of interest are studied. The main result lies in the fact that all 
their probability distributions obey a central limit theorem. 

1 Introduction 

In the literature coloured hard-dimers are applied in the framework of causally 
triangulated (2 -|- l)-dimensional quantum gravity. It was proved in [Ij, by 
using special triangulations of spacetime, that the generating function of the 
one step propagator depends on that one of coloured hard-dimers. For any 
configuration ^tvi of length A^, of blue and red sites, a coloured hard-dimer is 
a sequence of blue and red dimers, that satisfy the "hardness" condition, i.e. 
they can not intersect. A dimer is an edge connecting two nearest sites of the 
same colour. 
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In the present paper we consider the "averaged" generating function for 
coloured hard-dimers, that is, the mean of generating functions over ah the 
configurations ^at, with fixed. In [2] we have found an exphcit formula for 
it together with estimates from above and from below, that are both expo- 
nential, for large N. In the following we study the probability distribution 
associated with the averaged generating function. Then we analyze the prob- 
ability distributions corresponding to some random variable (r.v.) of interest. 
In particular the number of dimers and the total length of them. It turns out 
that the role played by the r.v. that measures the total number of dimers and 
single points (i.e. sites not occupied by dimers) is very important. We prove 
that the r.v. total length of dimers is binomial with parameters — 1 and 
|. Moreover we see that, even though the number of dimers has an unknown 
probability distribution, we are able to estimate its mean and variance asymp- 
totically, by using some recursive formulas that relate its moments with that 
ones resulting from the dimers' length. Although the dimers' number dis- 
tribution is not binomial, its variance is of order A^, as A^ — > oo, as in the 
binomial case. 

The main result of the present article is a local Central Limit Theorem 
(C.L.T.), for A^ large enough, for the joint probability distribution correspond- 
ing to the number of dimers and that one of dimers and single points. The 
limit distribution is a bivariate gaussian distribution with correlation coeffi- 
cient equal to — Hence a C.L.T. holds, as A^ ^ oo, also for the marginal 
probability distribution related to the dimers' number. 

The paper is organized as follows. In section 2 we define the probabil- 
ity distribution associated with coloured hard-dimers, through the averaged 

generating function, and find an exact expression for its normalizing constant 

, o> TV— 1 

Cat = (§) • Moreover we recognize the right probability distribution for 
the length of dimers. In section 3 we calculate the first two moments of 
the dimers' number. Finally in section 4 we prove a C.L.T. for the dimers' 
number. 

2 Coloured hard-dimers and probability 
distributions 

Given a sequence ^tv of length A^ of blue and red sites on the one-dimensional 
lattice 1i, one defines a dimer to be an edge connecting two nearest sites of 
the same colour, that characterizes the dimer colour. A sequence of coloured 
and non overlapping dimers in turn yields a "coloured hard-dimer". In Fig.l 
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an example of a coloured hard-dimer is given. 



Figure 1: A hard-dimer, = 13 



As described in [Ij and one introduces the generating function associ- 
ated to coloured hard-dimers on S,n 

Z^^{u,v,w) = ^n"''(^)t;"'-(^)ti;"'"-(^), u,v,w G (0,oo), 
D 

where D is a hard-dimer on ^n, nh{D) and nr{D) indicate the number of 
blue and red dimers respectively on D, and nbr{D) the total number of sites 
within each dimer (this means the sites having a colour different from the 
colour of dimers containing them). Moreover, let us define by 76 (-D) and 
7,.(L>) the number of blue and red sites of D respectively, not occupied by 
dimers ("single points"). Then the following constraint 

2nh{D) + 2nr{D) + ni^{D) + 7^(1?) + 7^(Z)) = N (2.1) 

holds. In the above example nf){D) = 1, nr{D) = 2, nbr{D) = 3, 7b(-D) = 3 
and 7r(-D) = 1. 

In [2] we studied the average of over all the configurations ^iv, i.e. 

Zn{u,v,w) = ^J2^iNiu,v,w) (2.2) 

and found estimates from above and from below for Z^v- In order to obtain 
them we proved an explicit formula for the mean Zjv, by using combinatorial 
tools 

TV [|] 

Zn{u,v,w) = 1 + 

t=l s=l 

([•] denotes the integer part). 

Now, let us define a family of probability spaces {Q,m,^n, PN)Nen>i- We 
choose r^Tv to be the set of all different hard-dimer configurations, where 




4 



by hard-dimer configuration we mean a particular sequence ^tv of length N 
together with a particular hard-dimer on ^jv- The u-algebra J^n is the set of all 
subsets of and for P/v we take the probability measure which is distributed 
uniformly on J^tv- Normalizing the function Zjy by ■^Ar|(tt,^,^)={i,i,i) then just 
gives the joint generating function of the random variables n^, n,., ribr, defined 
on (QjV) -^Afi -P/v)i which count, for each hard-dimer configuration, the number 
of blue, red dimers and the total number of sites within dimers, respectively. 

The random variable rib + n-r (characterized by the index s in formula 
(j2.3p ) gives the number of dimers for each hard-dimer configuration. Corre- 
spondingly, the random variable N — jb — Ir (indexed by t in ()2.3p ) denotes 
the number of sites occupied by dimers. The above formula (j2.3p is obtained 
by fixing first the number of blue dimers {rib), I'ed ones {rir) and single points 
(76 and 7r), i.e., n;, + + 76 + 7r, indexed by N — t + s, without assign- 
ing any length to dimers, even though the total length of dimers, defined as 
rib + rir + ribr, (indexed by t — s) is given because of (12. ip . This gives a factor 

(rib + rir + -fb + IrV- 
riblrirl 7fe!7r! 

Then, in the non-trivial case where rib + 7^ 0, one can assign a length 
to each dimer, taking into account that the total length of them is given by 
rib + rir + nbr, contributing another factor ("■ij+"'-+"i>'-~iy gy gumniing over 
7b, rib and then over s and t, we get the expression (12. Sp . See [2] for the details. 

In the present paper we want to analyze the random variables + 
and rib + rir + nbr, more precisely their probability distribution. We prove 
that the random variable rib + rir + nbr is binomial with parameters — 1 
and ^ and, hence, a Central Limit Theorem (C.L.T.) holds, for large (De 
Moivre-Laplace's Theorem). In the case of Ub + rij. we see that the probability 
distribution is unknown, but, fixing the random variable n;, + + 7b + 7r, 
which in turn is binomial with parameters — 1 and |, because of (12. ip . the 
conditional distribution of Ub + rir is hyper geometric. Moreover, for N large 
enough, a C.L.T. holds for the joint probability distribution of rib + nr and 
r^b + + 76 + 7r and, hence, also for the distribution of rib + ^r- For the proof 
of these results the random variable nb + n,. + 7b + 7r plays a very important 
part. 

Evaluating the averaged generating function at the point u = v = w = 
1 we derive the normalizing constant of the probability measure P/v associated 
to coloured hard-dimers 
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This also shows that the joint probabihty distribution related to the r.v.s 
Uh + n-r and 2(n{, + n^) + rihr, more precisely 

Pn{s, t) = Pn {rib + nr = s; 2{nb + rir) + ntr = t) (2.5) 

is given by 

The main result of this section is an explicit formula for the normalizing 
constant C^r, that holds for any N. It is obtained by using combinatoric 
arguments. 

Theorem 2.1 For any N , the following formula 

= (-) (2.7) 

holds. 

Proof: Consider the following change of variables t' = N — t,k = N — t + s, 
so that (j2.4p becomes (note that after changing the variables, we shall rename 
t' again by t) 



( [f ] k-l N-1 k-1 

Cn = 1 + 

\^k=l t=0 t=2k-N 



Note that in formula (|2.8|) only the combinatorial coefficients depend on t. 
Moreover, the combinatorial coefficients of the first sum in (j2.8p (^) {'^Zt-i) 
yield a non-normalized hyper geometric distribution with parameters — 1 
(population size), k (number of successes in the population) and k — 1 (sample 
size). Therefore summing over t we get 

[—1 f— 1 

^_ ^ /iV - - 1\ _ M ^_ /AT - 1\ 

2^2N-k2^(t)(k-t-l ) "^2^-n j ^ ' 

k=i t=o ^ ^ ^ ^ k=l ^ ^ 

Analogously, performing the change of variable t' = t — 2k + N, we get in the 
second sum of (j2.8p again a non-normalized hyper geometric distribution with 
parameters N — 1,N — k — 1,N — k — 1. Summing again over t we get 
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1^ 2N-k 1^ [t)\k-t-l 

fc=[f ]+l t=2k-N ^ ^ ^ 

^^^-V k \(N-k 
fc=[f ]+i t=o ^ ^ ^ 



E ^U-i) ^'-''^ 



with the convention that the binomial coefficient (^) = if fc > n. Putting 
together the last terms in (j2.9p and (j2.10p we get a binomial formula 

k=l ^ ' k=0 ^ ^ ^ 

□ 

An important Corollary of the previous Theorem is that the probability 
distribution of the r.v. nj, + + 76 + 7r (total number of dimers and sin- 
gle points) related to the index k in (j2.8p has a binomial distribution with 
parameters — 1 and |. 



Corollary 2.1 The probability distribution of the r.v. rih + rir + + 7r is a 

binomial distribution with parameters N — 1 and |. 

Proof: Performing the same change of variables as for (j2.8p and summing 
the probability distribution with respect to t, as in (|2.8p we get 

,2\^-ViV-l\ 1 

PNyrib + nr + ^h + lT = k) 



k-l 2^-^ 



N-l\ f2\^-^ fl^ 



k-lj \3J V3. 
for 1 < < - 1. For k = N it holds n?, + = 7f, + 7^ = A. Hence 

1 /2^^"^ 

PNirih + nr + 7f, + 7r = A) = PNiUb + nr = 0) = — = ( - 
by ([22]). □ 
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Corollary 2.2 For large N a C.L.T. (De Moivre- Laplace) for the r.v. Uf, + 

+ 76 + 7r holds 

^ _(fe-§(iv-i))' 
PivK + nr + ^b + lr = k) = . e t(^-i) (1 + rfc(iV)) 

where limjv^oo ^ki^^) = 0, uniformly with respect to k, for any k such that 



y = [k — ^{N — 1)] / Y |(-^ — 1) belongs to a finite interval {—A, A). 

Analogously one can find a similar result for the r.v. ni, + n.^ + n^r (total 
length of dimers). 

Corollary 2.3 The probability distribution of the r.v. rib + rir + f^br is a 
binomial distribution with parameters N — 1 and |. 



Proof: Taking into account (j2.ip one has 

PN{nb + nr + ribr = h) = PNim + rir + jb + 7r = N - h) 
N-l\ ^-^-^ 



h J \3J V3. 

for < /i < iV - 1. □ 

Corollary 2.4 For large N a C.L.T. (De Moivre- Laplace) for the r.v. rib + 

rir + Ubr holds 

_(h_l(iV-l))2 

PN{nb + nr + Ubr = h)= , e i^"-'^ (l + r'^{N)) 

where limjv_^oo '"^ (-^) = 0, uniformly with respect to h, for any h such that 



[h - \{N - l)]/v/|(A^ - 1) belongs to a finite interval {-B,B). 



3 Number of dimers: moments 

In the present section we investigate the distribution of the random variable 
Ub + Ur, more precisely, we calculate the first two moments of it, by using 
recursive asymptotic formulas that depend on the first two moments of the 
binomial r.v. nf, + + Ubr, studied in the previous section. Starting from 
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the averaged generating function Zjv, defined in (j2.3p . we rescale it by the 
normahzing constant Cat, calculated in the previous section 

Zn{u,v,w) 

Zn[u,v,w) = = 

Cat 

Therefore we calculate the moments of Uf, + through its derivatives (see 
[3]). By symmetry of the variables u and v in (j3.ip we have 

EN{nb) = -q;^Zn{u,v,w\u,v,w)={i,i,i) = 
d 

-^Zm{u,V,w)\(^u,v,w)=(1X'^) = EN{nr) = 

iv-i TV^^ g/2 /iV-t + g\ 

We indicate by the mean of random variables with respect to the proba- 
bility measure Pn- Therefore 

^ . ^ /2\^-iA^ s [N -t + s\ft-s-l\ 

iV [|] 

= EE«^^("'^) (3.2) 

t=l s=l 

We could deduce the above formula (j3.2p by noting that the index related to 
the r.v. rifo + is s. Analogously we find the formula for the mean of the r.v. 

Eninb + nr + Ubr) = EN{nb + n^) + EM{nbr) 

and 

dw 



d 

Ewiribr) = ■^Zn{u,v,w)\^u,v,w)={i,i,i) 



2X7V-1 t _ 2g /jV - t + g\ /t - g - 1 
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N I2I 
t=l s=l 

Hence 

i.„(n. + „. + n.) I + (3.3) 

In the next proposition we prove a recursive asymptotic formula, that hnks 
the mean of the r.v. + rir to that one of nh + nr + Uhr, for N large enough. 

Proposition 3.1 The following asymptotics 

N 2 
EN{nb + nr) X y - EN{nb + Ur + n^^) + - EN-i{nb + + Uhr) (3.4) 

holds, for N large enough. 

Remark 3.1 From Corollary 2.3 we have that under P^, n;, + + Uhr ~ 

B[N — 1, |) (binomial distribution), so that under Pn~i, n;, + + n^r ~ 
B{N-2,l). Hence 

, N-1 

ENynb + nr + Ubr) = — - — 
N -2 

EN~l{nb + Ur + Ubr) = — - — (3.5) 
From Remark 3.1 the next Corollary follows: 
Corollary 3.1 For large N , 

, 2N -1 , , 

ENint + nr) X (3.6) 



Proof: Applying the formula (j3.4p and taking into account (|3.5p we obtain 

, N -1 2 N -2 2N -1 

i?Mn, + n.)x--^ + -^ = ^— 

□ 
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Remark 3.2 By identity h2. 1\) and from hS. 5|) and 13. 6\) we are able to cal- 
culate asymptotically the single point number's mean. In fact 

27V - 1 - 1 4 , 
N = -(N + 1) 

Note that if we consider only the first order of the asymptotics with respect to 
N , we have 

4 

ENhb + 7r) ^-^^^ 2EN{nb + Ur) 

y 

that is, for the present model the expected number of single points is asymp- 
totically twice the expected number of dimers. Moreover, fixing the number 
of single points, the conditional probability distribution of jb (jr) is binomial 
and symmetric. 

Proof of Proposition 3.1: From (j3.2p we have 

(3.7) 

We applied in (j3.7p the identity 

Moreover, by Pascal's identity 

N -t + s-1^ _ ^ -t + s^ _ ^ -t + s-1 

the mean EN{nb + Uj-) becomes 



7V-1 N h] 



2*^ 

t=l s=l 



N -t + s\ /N -1-t + s 



't- s-V 
s - 1 

(3.8) 

Taking into account the definition of the probability distribution P/v in (|2.6p , 
the formula (13.81) becomes 



^7v(nfe+n,) = ^^(A^-i+s)P7v(s,t)-- ^ ^{N-t+s)PN-i{s,t) (3.9) 



3 

t=l s=l t=l s=l 
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In fact the normalizing constant for Pn-i is Cn^i = (|) 
The first sum in ()3.9p gives 



instead of Cat. 



N 



N 



^ E E ^^(«' - E E(* - ^)^^(^' 



t=l s=l 



t=l s=l 



N 



N-l 



EN{nb + nr + Ubr 



N - E^inb + nr + mr 

by (IMI). 

Analogously for the second sum in p.9p one has 



Af-l [|] N-l [|] 

^ E E^^-i("'*) - E T.(*-s)PN-i{s,t) 



t=l s=l 



t=l s=l 



(3,10) 



3 



N~2 



+ -EN^i{nb + Ur + ribr 



2 2 

-A^ + -EN-lim + nr + Ubr) 



(3.11) 



Both formulas (|3.10p and (|3.1ip give (j3.4p and Proposition 3.1 is so proved. 

□ 

In order to find the variance of Uf, + we need an analogous recursive 
formula for the second factorial moment of nft+n^ and that one of nb+rir+ribr, 
whose distribution is well known. By ()3.2p the second factorial moment of 
nb + rir + Ubr is of the form 



EN[{nb + nr){nb + Ur - 1)] 
TV [|] 

Y,Y. s{s-l)PN{s,t) 



(3.12) 



t=l s=l 

Analogously it easy to see that the second factorial moment of nb + rij. + Ubr 
is 

EN[{nb + nr + nbr){nb + Ur + Ubr - 1)] = 



N 



J2^{t-s){t-S-l)pN{s,t) 



(3.13) 



t=l s=l 
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In fact, the indices related to the r.v.s rif, + and n?, + + n^j. are s and 
t — s respectively. 

We generalize the asymptotic recursive formula for the first moments, given 
in Proposition 3.1 to the second factorial moment of Uf, + rir in terms of the 
first two factorial moments of + + Ubr- 

Proposition 3.2 For N large enough, the following asymptotic recursive for- 
mula 

M N(N-l) 

EN[{nb + nr){nb + - IjJ x h 

y 

4 4 
-2{N-l)[EN{nb+nr+nbr)-- EN-i{nb+nr+nbr) + - EN-2{nb+nr-\-nbr)] + 

4 

EN[{nb+nr+nbr){nb+nr-\-nbr-l)]--EN^i[{nb-\-nr+nbr){nb+nr-\-nbr-l)] + 
4 

-EN_2[{nb + nr + nbr){nb + + Ubr - 1)] (3-14) 

y 

holds. 



Proof of Proposition 3.2: From (j3.12p and taking into account the proof 
of Proposition 3.1, we have 

EN[{nb + nr){nb + Ur - 1)] = 

2\^~^ {N - t + s){N - t + s - I) fN -t + s-2\ ft-s-l 

3/ 2*^^ V •s-2 )\ s-l 

^ t=l s=l ^ ^ ^ 

because of the identity 

^^fN-t + s\ ,,,, ^.fN-t + s-2 

s{s-l)i^ ^ j ={N-t + s){N-t + s-l)i^ 

Considering the expressions ()3.12p and (j3.13p of the second factorial moments 
of Ub + Ur and nb + nr + Ubr respectively, we rewrite (A'' — t + s){N — t + s — l) 
in a suitable form 

{N-t+s){N-t + s-l) = N{N-l)-2{N-l){t-s) + {t-s){t-s-l) (3.15) 

Moreover, as in Proposition 3.1, we apply (this time twice) Pascal's formula 

N-t + s-2\_fN-t + s-l\ /N -t + s-2 
s-2 ) ~ [ s-l j~\ s-l 
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N-t + s\ ^/N-t + s-l\/N-t + s-2 



(3.16) 



By (|3.15p and (j3.16p the second factorial moment of Uf, + rir becomes 

ENiirib + nr){nb + - 1)] = 

7V-1 N [|] 



^ Y^[N{N - 1) - 2{N -l){t-s) + {t-s){t-s- 1)] X 

t-s-l\ 1 



t=l s=l 



N-t + s\ ^/N-t + s-l\^/N-t + s-2 



N 



1 /2* 



^ Y^[N{N - 1) - 2{N - l)(t -s) + {t- s){t - s - l)]PN{s,t)+ 

t=l s=l 

-3 E - 1) - - - + (t - - S - l)]PjV-l(5, t)+ 

t=l s=l 

iV-2 [|] 

+ 9 E E[^(^ - 1) - - l)(t - + (t - - s - l)]Pjv-2(5, t) (3.17) 



t=l s=l 



where the normalizing constant for Pn-2 is Cn-2 = (|)^ ^• 
The first sum of the right hand side of (j3.17p becomes 

N [|] N [|] 

N{N -i)Y.Yl ^^(«' - - 1) E E(* - *)^^(^' 



t=l s=l 
N [|] 



t=l s=l 



Y^Y,{t-s){t-S-l)PN{s,t) 



t=l s=l 



N{N - 1) 



1 



N-l 



2{N - l)EN{nb + nr + nbr) + 



EN[{nb + nr + nbr){nb + nr + nbr-l)] N{N -1)-2{N -l)EN{nb + nr + nbr) + 
EN[{nb + nr+ nbr){nb + Ur + Ubr - 1)] (3.18) 

As above the second sum in (j3.17p is of the form 

N-2' 



--N{N-l) 



1 



+ -[N - l)EN-i{nb + nr + nbr) + 



14 



--EN-i[{nb + nr + nhr){ni, + rir + Ubr - 1)] >i 

4 8 
--iV(iV - 1) + -(iV - l)EN-i{nb + nr + nbr)+ 

4 

- -EN-l[{nb + nr + nbr){nb + Ur + Ubr - 1)] 



Finally the last sum in (j3.17|) gives the following contribution 
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A^(A^ - 1) 



1 



(3.19) 



-{N - l)EN~2{nb + nr + nbr) + 



+-EN-2[{nb + nr + nbr){nb + rir + Ubr - 1)] 
4 8 

-N{N - 1) - -{N - l)EN-2{nb + nr + nbr)+ 
4 

-EN_2[{nb + nr + nbr){nb + + Ubr - 1)] 



(3.20) 



Putting together p.l8p - (|3.20p we get (j3.14p and Proposition 3.2 is so proved. 

□ 

Now we are able to calculate the second factorial moment and the variance 
of Ub + rir 

Corollary 3.2 For large N , the following asymptotics 



EN[inb + nr){nb + - 1)] x — (iV - l)(iV - 3) (3.21) 

oi 



holds. 



Proof: Since rib + '^■r + ^frr ~ B{N — 1, i) for sequences ^at of length A^, one 
can easily find its second factorial moment 



EN[{nb + rir + ribr){nb + rir + n-br - 1)] 
Analogously we have 

EN-iiirib + nr + nbr){rib + + ribr - 1)] 
EN_2[{nb + rir + ribr){rib + Ur + Ubr - 1)] = 



(A^-l)(iV-2) 
9 



[N -2){N -?>) 

' 9 
(iV-3)(iV-4) 

9 



(3.22) 



(3.23) 
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From (j3.5p . (j3.22p . (j3.23p and Proposition 3.2 we obtain the second factorial 
moment of Uh + rir 



N(N-l) 

^ '--2iN-l) 



EN[{nb + nr){nb + Ur - 1)] x 

N-1 4.N-2 AN 



3 3 3 9 



+ 



{N-l){N-2) 4{N-2){N-3) A{N-3){N-4)_ 
9 3 9 9 9 ~ 

_ 2 _ ^ _ _ 2) _ ^(^ _ _ 3)+ 

-^{N - 2){N - 3) + l(iV - 3)(iV - 4) = 

/V — 1 32 

__[Ar + 9(iV - 2) - 6{N - 1)] - -(iV _ i)(N - 3) = 

l(Ar_i)(iv_3) 



□ 



Corollary 3.3 For N large enough 



VarN{nb + nr)^^{3N + l) (3.24) 

The symbol Var^ indicates the variance associated to probability measure P^r. 
Proof: One has that 

VarN{nb + nr) = EN[{nb + nr){nb + nr - 1)]+ EN[nb + nr] - (i?Ar[nfe + n^])^ x 

— [4(iV - 1){N - 3) + 9(2iV - 1) - {2N - 1)^] = 
81 

|f(3iV + l) 

□ 

Remark 3.3 Note that in the variance formula {3.24^ for Ub + the second 
order term with respect to N disappears, so that the accuracy of the first order 
terms with respect to N is important. Nevertheless, from the asymptotics for 
the mean and the variance ofub+n^ ( liS. 6\) and \3.24-^ ) one can deduce that the 
distribution of Ub + is asymptotically not binomial. In fact EN{nb + n^) ^ 
^N and Var]\f{nb + rir) §fN ^ | • \N . In the next section we prove that it 
is asymptotically gaussian, for large N, i.e. a C.L.T. holds. 
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4 Central Limit Theorem for dimers' num- 
ber 

In the present section we study the asymptotic distribution of the dimers' 

number, in particular we prove a C.L.T. for the joint probability distribution 

of the total number of dimers and single points (n^ + + 76 + 7r), analyzed 

in Section 2, and the number of dimers (nj, + n^). The limit distribution is 

a bivariate gaussian distribution with correlation coefficients equal to — 7=. 

v3 

The proof is a generalization of De Moivre-Laplace's Theorem. 
Theorem 4.1 A C.L.T. holds for the joint probability distribution 
PN{nb + Ur = s;nb + rir + + 7r = k) = 

1 _3(,2+2^,j,+j,2)^^^^^^^^^^^^ ^^^^^ 



(2vr)2 (Aat) (Iat) (I) 



where 



y = -7=r- (4-2) 



9 3 

Moreover linijy^oo = 0, uniformly with respect to x and y, defined in 
14.^ , belonging to finite intervals {—C,C) and {—A, A) respectively, with C 
and A positive real constants. 

Proof: As in Section 2 we perform the change of variables k = N — t + s and 
s = s on Piy{s,t), defined in (j2.6p . The index k indicates the total number of 
dimers and single points. Then the probability PN{s,t) becomes 



PN{nb + Ur = s;nb + Ur + jb + Ir = k) = Pn{s, N - k + s) = 

(4.3) 



J \ s — 1 / V^y 

Taking into account that the indices s and k are both of order N, we can 
forget in (14.30 the constants, i.e. k — 1 ^ k and N — k — 1 x N — k, as 
iV ^ 00. 

As in De Moivre-Laplace's Theorem, we apply Stirling's formula to the 
binomial coefficients in (j4.3p . In the present model we have two binomial 
coefficients instead of one, so that the calculus becomes heavier than in De 
Moivre-Laplace's Theorem. We write 

k\ /N-k\ f2\^ 



s I \ s J \ 3 J \3 
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k] {N-k)l f2Y fl 



/I \ N-k 



s\{k - sy. sl{N - k - s)] \3j 



s^{k - s)*=-^ s^(iV - k- s) 

p^N,3,k 



N-k-s 



where 
and 



r, s(k—s) /o s(N—k—s) 
^'^^^V^^ N-k 



^N,s,k = Afc — As — Xk-s + AtV-A: — — AtV-A; 



< An < T7— , for any n G N 



(4.4) 



12n + 1 - 12n 
Then we rewrite (14.411 as 



3s J \3{k-s)J \ 3s J \3{N-k 



(4.5) 



In order to find an asymptotics for ()4.5p we recall that the r.v. nft + nr+7fe+7r 
has binomial distribution B{N — 1, |) (Corollary 2.1) and hence Ejy[{nh + nr + 
lb + lr)-l{N -l)^ In and VarN{ni, + + 7,, + 7^) x |(iV - 1) x . 
Moreover in the previous section we proved that En (rib + n^.) x ^^q~^ x 
and that V^arAr(n6 + ra^) x ^(3iV + 1) x ^A^. 

According to De Moivre-Laplace's Theorem we normalize the variables 
nb + Tir + lb + 7r and nb + Ur, by using the asymptotics of the moments and 
performing then the change of variables (|4.2p so that 



s = |iV + for X G (-C, C) 



k = \N + ^y ioiye{-A,A) 



(4.6) 



Taking into account (j4.6p . we consider the first factor ^at^s^a; in ()4.5p with 
respect to the variables a; and y, in particular 



XN,s,k < 
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12 f |iV + Vmy\ 12UN + ^x) +1 12 f ^iV - + ^y) + 1 



^12(f-^,) 12(f-^x-4^,)+l ^'-'^ 

Since x, y belong to bounded intervals one can estimate XN,s,k uniformly from 
above with respect to x and y and we get a bound of order that doesn't 
depend on x and y. Analogously, we can find a uniform lower bound of order 
jj. Therefore 

with r2']j{N) 0, as N ^ 00, uniformly with respect to x and y. 
The second factor in (j4.5p is Biy^s,k 



B 



N,s,k 



From (j4.6p we estimate 

s{k-s) + {^^-^- + ^y, 



k {iN + ^y 

M_l_ /-6.2IIM_ /^2;_l3 /J_„, 
4 ^ \l N 2 ) Y \J N i\J N y 

\ ^ V ^ 2 

^7v(l + rg)(iV)j (4.8) 

with rx,y{N) ^ 0, as — > 00, uniformly with respect to x and y, since 
x G (— C, C) and y G (— ^, so that the upper and lower bounds of B]\[^s,k 
do not depend on x and y. 
Analogously 

s{N -k-s) 2 



N -k 27 

,(3) 



A(l + rg)(A)) 



with rx,y{N) ^ 0, as A^ ^ 00, uniformly with respect to x and y. Therefore 

BN,s,k r- (4-9) 



9 3 V 3 
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as ^ oo. 
Note that 



B 



N,s,k 



2TTy/VarN{nb + nr)^^VarN{nb + n,. + 7;, + 7r)\/l - 



where p = is the correlation coefficient, whose sign will be determined 
later. 

Finally we consider the logarithm of the last factor CM,s,k in (|4.5p 



lnCAr,s,fc = -sin 
3s 



(fc-s)ln('^^^'l + 



-s In 



2{N - k) 



{N -k-s)\ii 



2k 

3{N -k-s) 



N -k 



i=l 



N,s,k 



We express now each term of the sum in ()4.10p ^, 
of X and y, defined in (14. 2p . We start with 



(4.10) 

i = 1,2,3,4 in terms 



N,s,k 



9 9 



6N 



-X In 



'2N 



^ + \/3 2; I In 1 + 



^/3 



2; - 2/ 



/2iV + y 



Since the last logarithm above is of the form ln(l + z), with z ^ 0, we can 
expand it around 2 = 0, ln(l + z) = z — ^ + 0(2;^), as z ^ 0. The same is 
true for each logarithm function present in any ^ ^, i = 1, 2, 3, 4. So C\j ^ ^ 
becomes, as ^ 00, 



C 



1 

N,s,k 



/2N 



^ + x){^/^ X - y] ( 2\/2iV - ^/3 x + 3y 



/2N + y 



(4.11) 



Analogously for ^ ^ 



CAf,s,fc = -(A: - s)ln 



3(fc - s) 
2fc 



20 



6N 



-x + 



/2N 



-y In 



N 



x + V2N y 



/2N 



(2V2N - X + 3j/) In [1 + 



-V3x + 



+ 



9 ^ ■ " ■ - - ■ - y- ■ 2{V2N + y)^ 
/2N (2V2N -^/Sx + Sy^ (VSx-y) (aV2N + x + 3y 



72 



^+yy 



Then C%^^^, 



CN,s,k = — sin 



3s 



2{N-k)_ 
( 2N + V6NX 



2 V6iV \ ^ 

-AT H —X In 

9 9 / \2N-2V2Ny 



/2N 



'2N + 



VS In ^1 + 



^/Zx + 2y 



'2N -2yj 
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(y2N-2yy 



It remains to see ^ ^ 



N VGN V2N \ 

:^ — X — y m 

9 9 3 ^1 

Vm ( [n 



N-k 

3 3 y ' 



(^yf -V3x-3yj Inh 



V3a; + 2y 



(\/f - \/3 ,r - 3?/^ (\/3 ,r + 2?/) + \/3 ,r) 

(V2iV-2y)' 



(4.12) 



(4.13) 



(4.14) 
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Summing the last term in (j4.1ip {Cj^^ ^) with that one in (|4.12p (C^^ we 
obtain 



C^N,s,k + C^N,s,k 



/2N 



2V2N -V3x + 3y){V3x-y 



(4.15) 



Summing the last term in (j4.13p (C^ ^ ^) with that one in (j4.14p (C^ ^ ^) we 
obtain 



CN,s,k + 



f2N (\/2iV + \/3x)(^/3x + 2yy 



6 



'2N - 2y 

Finally, the main contribution of (j4.15p and (|4.16p is 



1=1 



N,s,k 



1 

12 



3 
4 



\/3 X - V 2 f \/3 X + 2y^ ^ 



2 r, 2 



(4.16) 



(4.17) 



Remark 4.1 iVoie i/iai the last term in j^^. j?] ) z/ 0/ i/ie /orm 

1 



2(1 



2pxy + y ) 



with /> = — -1=. IFe yet so a bivariate gaussian distribution with correlation 

coefficient equal to *-e., the r.v. rib + and + + 76 + 1r are 

negatively correlated. 

From the previous theorem we can deduce the following result. 
Corollary 4.1 For large N, a C.L.T. for the r.v. Ub + holds, i.e. 

1 



P{nb + nr = s) 
for any s defined in (|^.g[). 



-e 27" 



(4.18) 
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Proof: From (j4.3p we have 

»'^|e)r:-r')(l)'"'(r"'^ 

/ iN-^A IN+^A iN-^x \ 

E + E + E Uiv(s,iV-A; + s) (4.19) 

^fc=f7V+4Hx fc=§iV-v|lA k=^N+^Aj 

In (|4l9]) . let us denote by A^), S2(A, N) and S3(A, iV) the first, second 

and third sum, respectively. Using Theorem 4.1 and the integral C.L.T. of De 
Moivre-Laplace [1| one finds for every e > an A^(e) such that 

\j:i{A,N)-Ii{A)\<^, i = 1,2,3 (4.20) 

for any N > N{e). The terms Ii{A) are given by 

1 _Ei 1 f -vl , 
-.e 2 / e 2 



27" 

-A r rA 



and the boundaries of the integrals f. are fixed as = J_^, = and 
= J^, respectively Note that K20h proves the statement. □ 
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